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SUMMARY 

Four formulas, for the nonzero altitude transformation from geodetic coordinates 
(geodetic latitude and altitude) to geocentric coordinates (geocentric latitude and geocen- 
tric distance) and vice versa, are derived. The set of four formulas is expressed in each 
of three useful forms: series expansion in powers of the earth's flattening; series expan- 
sion in powers of the earth's eccentricity; and Fourier series expansion in terms of the 
geodetic latitude or the geocentric latitude. The error incurred in these series expan- 
sions is of the order of one part in 3 x 10^. 

INTRODUCTION 

In astrogeodetic work one frequently needs to transform from geodetic coordinates 
(geodetic latitude and altitude) to geocentric coordinates (geocentric latitude and geocen- 
tric distance). The inverse transformation is also frequently needed. 

Before the days of air and space travel, locations of interest remained on the sur- 
face of the earth. Because the altitude above mean sea level of locations on the earth's 
surface is much smaller than the earth's radius, a number of transformation formulas 
were derived for zero altitude (above mean sea level) and for the geocentric distance 
equaling the earth's radius. A formula for the geocentric latitude as a function of the 
geodetic latitude for zero altitude is given as a Fourier series expansion in references 1 
to 5. A formula for the geocentric distance as a function of the geodetic latitude for zero 
altitude is given as a Fourier series expansion in references 3 to 6 and as an analytical 
expression in references 1 to 3. Also, a formula for the geodetic latitude as a function 
of the geocentric latitude for the geocentric distance equaling the earth's radius is pre- 
sented as a Fourier series expansion in references 4 to 6. 

After the arrival of the days of air and space travel, locations of interest shifted to 
above the surface of the earth. A number of transformation formulas have been derived 
for any altitude and geocentric distance. Analytical expressions for the geocentric 



distance as a function of the geodetic latitude and the altitude are given in references 7 
and 8. A formula for the geodetic latitude as a function of the geocentric latitude and the 
geocentric distance is given as a Fourier series expansion in reference 9 and as an 
approximate expression in reference 10. Several analytical expressions (three in ref. 9, 
one in ref. 10, and one in ref. 11) are presented for the altitude as a function of the geo- 
centric latitude, the geocentric distance, and the geodetic latitude. An approximate 
expression for thp altitude as a function of the geocentric distance and the geodetic lati- 
tude is given in reference 6. All six of these formulas, however, possess the same 
shortcoming; that is, for altitudes close to the earth's surface, the altitude is expressed 
as the difference of two comparable numbers and, consequently, the accuracy is impaired. 

In present-day astrogeodetic work, one needs a complete set of transformation 
formulas whose accuracy remains consistent for all altitudes. In this paper, formulas 
for the geodetic latitude and the altitude as a function of the geocentric latitude and the 
geocentric distance and vice versa are derived for nonzero altitudes. For the reader's 
convenience, the four formulas are cast as series expansions in powers of the earth's 
flattening, as series expansions in powers of the earth's eccentricity, and as Fourier 
series expansions in terms of the geodetic latitude or the geocentric latitude. 

Some astronomical applications (such as the laser ranging work in refs. 12 and 13) 
require the extra precision in coordinate determination that is possible with satellite 
geodesy techniques. However, for most environmental and space sciences applications 
(such as the tracking of space vehicles, satellites, and missiles in refs. 6, 10, and 11 and 
the magnetospheric physics work in refs. 7, 8, and 14), the formulas derived in this paper 
should be sufficiently accurate and very useful. If improved accuracy is required, addi- 
tional terms may be derived by use of the procedures described herein. 

SYMBOLS 

a semimajor axis of ellipse or equatorial radius of planet, arbitrary units 

a^ag^i'^' expansion coefficients for e in powers of f 

b semiminor axis of ellipse or polar radius of planet, arbitrary units 

bpbg,^'^' expansion coefficients for difference between h + 1 and p in 

powers of f 

C distance from point on ellipse, along perpendicular to ellipse, to intersection 

point with minor axis, units of a 
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e eccentricity 

f flattening 

h altitude, units of a 

S distance from point on ellipse, along perpendicular to ellipse, to intersection 

point with major axis, units of a 

x,y coordinates of point on ellipse, units of a 

€ difference between geocentric latitude and geodetic latitude 

<p geodetic latitude 

c p 1 geocentric latitude 

p geocentric distance, units of a 

ANALYTICAL FORMULATION 


Relations Pertaining to an Ellipse or to a Planet 
The equation of an ellipse is 


x 2 

a 2 



( 1 ) 


where a is the semimajor axis (or the equatorial radius of the planet) and b is the 
semiminor axis (or the polar radius of the planet). The polar radius b is related to 
the equatorial radius a by 

b 2 = a 2 (l - e 2 ) (2) 


The quantity e is the eccentricity, which is represented by the fraction of the major axis 
that lies between the two foci of the ellipse. 

The flattening f, which is the compression of the planet’s spheroid from an exact 
sphere, is defined by 

f _ a - b 
* ~ a 
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or 


b = a(l - f) (3) 

Comparing equations (2) and (3) yields the following equation relating the eccentricity e 
to the flattening f: 

1 - e 2 = tl - f) 2 (4) 

Substituting equation (3) into equation (1) results in 


y 2 . 

a 2 a 2 (l - f)^ 


(5) 


Differentiating y in equation (5) with respect to x renders the slope dy/dx at the 
point x,y on an ellipse. This slope is expressed as 


&=-(l -f) 2 ^ 
dx y 


( 6 ) 


Now, 


^ tan <p= -1 (7) 

dx 

where the angle <p is the geodetic latitude, shown in figure 1, which is the angle which 
the perpendicular from the point in space of interest to the earth's spheroid makes with 
the earth's equatorial plane. The earth's equatorial plane is the plane through the great 
circle that is halfway between the earth's North and South Poles. Substituting equation (6) 
into equation (7) furnishes 


- = (1 - f) 2 tan (j) 
x 


Substituting equation (8) for y into equation (5) produces 


x 

a 


1 + (1 - f) 2 tan 2 $| 


— i 
0 


= 1 


(8) 


4 



or 


a cos 0 

jcos^ + (1 - f )2 sin^^j^ 2 

Substituting equation (9) into equation (8) leads to 


(9) 


y _ a(l - f) 2 sin 0 

[cos + (1 - f)^ sin^^j ^ 

From figure 1 one sees that 

x = C cos 0 

y = S sin 0 


( 10 ) 


(ID 

(12) 


where C and S are the distances from the point x,y along the perpendicular to the 
ellipse to the intersection points with the minor and major axes, respectively. 



Figure 1.- The geodetic coordinates 0 and h and the geocentric 

coordinates 0’ and p. 
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Comparing equations (9) and (11) establishes that 


C = 


|cos 2 0 + (1 - f) 2 sin 2 0) 


1/2 


( 13 ) 


Comparing equations (10) and (12) gives 

« a(l - f) 2 

{cos 2 0 + (1 - f) 2 sin 2 *^)^ 

Comparing this equation with equation (13) yields 

S = (1 - f) 2 C (14) 

Expanding equation (13) in a binomial series expansion, setting a equal to unity, 
and omitting the terms of order f3 and higher result in 

C = 1 + 4f + 7 wf 2 -i(f+ f 2 )cos 20 + ^ f 2 cos 40 (15) 

2 16 2 16 

Note that, because of the small value of the earth's flattening f, the error incurred 
because of the omission of the terms of order f 2 and higher is of the order of one 
part in 3 x 1C)7. 

The use of equation (14) for S renders 
C - S = C - (1 - f) 2 C 
= (2f - f 2 )C 

3 

Substituting equation (15) into this equation and omitting the terms of order f and 
higher yield 


C - S = 2f - f 2 cos 20 (16) 

From figure 1 one observes that for an arbitrary point in space 

p cos 0' = (C + h)cos 0 (17) 
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p sin 0* = (S + h)sin 0 


(18) 


The distance h is the altitude, which is the perpendicular distance from the point in 
space of interest to the earth's spheroid. The distance p is the geocentric distance, 
which is the distance from the point in space of interest to the earth's center; and the 
angle 0' is the geocentric latitude, which is the angle that this radial line from the 
earth's center makes with the earth's equatorial plane. 

The difference between the geocentric latitude 0' and the geodetic latitude 0 
is denoted by e and is shown in figure 1 . This difference is expressed as 

e = 0 - 0' (19) 

Since the geodetic latitude differs by a small amount from the geocentric latitude because 
of the earth's flattening f, which is a small value, then 


sin e = e 


The quantity e is of the same order of magnitude as the earth's flattening f. The 
following Taylor series expansions where the terms of order and higher have been 
omitted will be used in the subsequent derivations of the transformation formulas: 


cos (0’ +e) = cos 0' - e sin 0’ - — cos 0’ 

Ct 


cos (0 - t) = cos 0 + e sin 0 - cos 0 


sin (0' + e) = sin 0' + e cos 0' - — sin 0' 

2 


cos(20' + 2e) = cos 20' - 2e sin 20' - 2e^ cos 20' 


sin(20’ + 2e) = sin 20' + 2e cos 20' - 2e^ sin 20' 


cos(4 0' + 4e) = cos 40' - 4e sin 40' - 8e^ cos 40' 
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Transformation Formulas in Forms of Series Expansions 
in Powers of Earth's Flattening 

Geodetic latitude .- Multiplying equation (17) by sin 0 and equation (18) by cos 0 
and then subtracting the second of these two resulting equations from the first produce 


p sin (0 - 0') = 4 (C - S)sin 20 (27) 

Ci 

The (^-dependence in equation (27) can be eliminated through the use of equation (19) to 
obtain 


p sin e = i (C - S)sin(20' + 2e) (28) 

Ci 

Substituting equation (20) on the left-hand side and substituting equation (25) on the right- 
hand side of equation (28) lead to 

pe = ^ (C - S)(sin 20' + 2e cos 2 0' - 2e^ sin 2$') (29) 

Using equation (19) to eliminate the 0-dependence from equation (16) and then using 
equation (24) establish that 

C - S = 2f - f 2 (cos 20' - 2e sin 20' - 2e 2 cos 20') 

Recalling that e is of the same order of magnitude as f and then omitting the terms 
of order f 2 and higher render 

C - S = 2f - f 2 cos 20' (30) 

Substituting equation (30) into equation (29) gives 

pe = I (2f - f 2 cos 20' )(sin 20' + 2e cos 20’ - 2e 2 sin 20’ ) (31) 

Since the quantity e is of the same order of magnitude as the earth's flattening f, 
one can assume that e can be expanded in powers of f. This power series expansion 
can be written as 

e = ajf + a 2 f^ (32) 
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where the terms of order f** and higher have been omitted and the coefficients a^ 
and a .2 are to be determined. 

Substituting equation (32) into equation (31) and then equating like powers of f 
yield the coefficients a^ and ag as follows: 

( 33 ) 


_ sin 20* 
1 P 


a 2 " 



sin 40* 


Substituting equations (33) and (34) into equation (32) yields 


(34) 



(35) 


Substituting equation (19) into the left-hand side of equation (35) proves that 


0 = 0 ’ 



(36) 


Therefore, equation (36) is the transformation formula, in the form of a series expansion 
in powers of the earth's flattening, for the geodetic latitude as a function of the geocentric 
latitude and the geocentric distance. 

Altitude .- Multiplying equation (17) by sin 0' and equation (18) by cos 0' and 
then subtracting the second of these two resulting equations from the first result in 


0 = C cos 0 sin 0* - S sin 0 cos 0' - h sin(0 - 0') 


or 


(h + l)sin(0 - 0') = (C - S)sin 0 cos 0' - (C - l)sin(0 - 0') (37) 


Subtracting equation (37) from equation (27) renders 

fp - (h + ljlsin(0 - 0') = (C - S)sin <#cos 0 - cos 0') + (C - l)sin(0 - 0') 

(38) 
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The $ -dependence in equation (38) can be eliminated through the use of equation (19) to 
obtain 

(p - (h + l)]sin e = (C - S)sin(0* + ejjcos^' + e) - cos 0|] + (C - l)sin e (39) 

Substituting equations (20), (23), and (21) into equation (39) furnishes 

(p - (h + l))e = (C - S)/sin 0' + e cos 0' - sin <p*\ e sin 0* - cos 0’^ + (C - l)e 


or 


p - (h + 1) = (C - S)^sin 0’ + e cos 0* - -y sin 0'^ ^-sin 0* - | cos 0*j + (C - 1) 
Omitting the e 2 term in this equation gives 

p - (h + 1) = -(C - SjJsin 2 ^ + sin 20* + ^-(cos 2 $' - sin 2 0’)J + (C - 1) 

Using equation (35) in this equation and omitting the terms of order f 2 and higher 
produce 


p - (h + 1) = -(C - S) jsinfy +|||- in p 2 ^ , ) f + - ^jsin 40'| 

|sin^20 j f2 ^ cos 20t _ S i n 20t^\ + (C - 1) 


i* > sin 2 0' 


1 

+ 2 


= -(C - S)< ( 1 ~ cos 20') + 3{ (1 - cos 40') 

2 8 p 


+ f* 


|(^ " ^ sin 4 0’ sin 20’ + 1 ( — p 2 ^ -) (cos 2 0’ - sin 2 <^') 


+ (C - 1) 


(40) 


The $- dependence in equation (15) can be eliminated through the use of equation (19) to 
obtain 


C-l = if + ^f^- i(f + f^)cos(20' + 2e) + ^ f 2 cos(40* + 4e) 
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Using equations (24) and (26) in this equation leads to 


C - 1 = I f + -A f 2 _ l(f + f 2 )(cos 20' - 2e sin 20' - 2e 2 cos 20’) 

2 16 2 

+ A f 2 (cos 4 0’ - 4e sin 4 0' - 8e 2 cos 40* ) 

16 

o 

Using equation (35) in this equation and omitting terms of order f and higher 
establish that 

c - 1 =| f +A f2 -I(f+ f 2 )[cos 2 0’ - 2f( s - A - 2 A jsin 2,] + A f 2 C os 40’ 

= I f + f2 - |( f + f 2 )[cos 2<t>' - f + A f 2 cos 40’ (41) 

Substituting equations (30) and (41) into equation (40) and omitting terms of order f 2 
and higher give 


p - (h + 1) = -(2f - f 2 cos 20 ’)|AjL-£^_ 2 £I) + 3f (Lz_£ 2§-MI)J 

+ | f + f2 - |(f + f 2 )[cos 20’ - f (1 ~ c ° s + A f 2 cos 40* (42) 

Just as one assumed that the difference e between the geocentric latitude 0' 
and the geodetic latitude 0 could be expanded in powers of the earth’s flattening f, one 
can assume that the difference between the quantity h + 1 and the geocentric distance p 
can also be so expanded. (Recall that, in this paper, the altitude h is in units of the 
equatorial radius a and a has been set equal to unity; hence, the quantity h + 1 
actually represents the quantity h + a.) The power series expansion for the difference 
p - (h + 1) can be written as 

p - (h + 1) = bjf + b 2 f 2 (43) 

where the terms of order f 2 and higher have been omitted and the coefficients bj 
and b 2 are to be determined. 
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Substituting equation (43) into the left-hand side of equation (42) and then equating 
like powers of f yield the coefficients bj and b£ as follows: 


, cos 2d>* - 1 
b l = 2 


■li-H) 


(cos 40* - 1) 


(44) 

(45) 


Substituting equations (44) and (45) into equation (43) proves that 


f 2 (46) 

Therefore, equation (46) is the transformation formula, in the form of a series expansion 
in powers of the earth’s flattening, for the altitude as a function of the geocentric latitude 
and the geocentric distance. 

Geocentric latitude .- The ^'-dependence in equation (37) can be eliminated through 
the use of equation (19) to obtain 

(h + l)sin e = (C - S)sin 0 cos (<£ - e) - (C - l)sin e (47) 

(Note that, just as the quantity h + 1 actually represents the quantity h + a, the quan- 
tity C - 1 actually represents the quantity C - a.) 

Substituting equations (20) and (22) into equation (47) yields 


h = (p - 1) + (' - e ° s 2 »') f + p - i)(l - COS m 


(h + l)e = (C - S)sin $^cos + e sin $ ~ ^ cos ~ ^ ~ 

= (C - S)6i2^4. + £ . $1 cos *1- (C - 1)< 


Using equations (16) and (15) in this equation results in 


(h + l)e = (2f - f2 cos 20)[^LM + e cos 2£) . cos 

' [l f + ll f2 ' H f + f2 ) cos 2ct) + n t2 cos 40 ] e 


(48) 
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Again, one can assume that the difference e between the geocentric latitude 0' 
and the geodetic latitude 0 can be expanded in powers of the earth's flattening f. This 
power series expansion can be written as 

e = aj' f + a 2 ' f 2 (49) 

where the terms of order f 2 and higher have been omitted and the coefficients aj ' 
and a 2 ' are to be determined. 

Substituting equation (49) into equation (48) and then equating like powers of f 
yield the coefficients a^' and a 2 ' as follows: 


, sin 20 
h + 1 


(50) 


a 


2 


sin 2<fr 
2(h + l) 2 


1. + — -]si 

(4(h + l) 2 4(h + l)J 


sin 40 


(51) 


Substituting equations (50) and (51) into equation (49) renders 

e „ fenMf + f sin 2» . r . l + 4 A f 2 

\ 1> + 1 / V 2(h + l) 2 (4(h * l) 2 4(h + 1)J j 

Substituting equation (19) into the left-hand side of equation (52) proves that 


(52) 


0 * = 0 + 


-sin 20\ . J -sin 20 f 1 

> h + 1 / \ 2(h + l) 2 [4(h + l) 2 



(53) 


Therefore, equation (53) is the transformation formula, in the form of a series ejq)ansion 
in powers of the earth's flattening, for the geocentric latitude as a function of the geodetic 
latitude and the altitude. 

Geocentric distance .- The 0'-dependence in equation (38) can be eliminated through 
the use of equation (19) to obtain 


Q) - (h + l)jsin e = (C - S)sin 00os 0 - cos (0 - e)J + (C - l)sin e (54) 
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Substituting equations (20) and (22) into equation (54) furnishes 

Q) - (h + 1 J) € = (C - S)sin 0 ^-e sin 0 + ^ cos <pj + (C - l)e 


p - (h + 1) = (C - S)sin <(> (-sin 0 + ^ cos <p) + { C - 1) 


* / 


= -(C - S)^sin 2 0 - c + (C - 1) 

The use of equation (52) in this equation produces 

p - (h + 1) = -(C - S) jsin 2 0 - 

+ I + lJ £ 

\2{h + l) 2 [4(h + l) 2 4(h + 1)J 


i-^jsin 4^f 2 j^^U(C-l) 


_ am 1 - cos 20) t (1 - cos 40) 

_V 'V 2 ‘ 8(h + 1) 


2(h + l) 2 


f- — k + - — ^jsin 40 

[_4(h + l) 2 4(h + 1)J 


+ (C - 1) 


Using equations (16) and (15) in this equation and omitting terms of order f 2 and higher 


lead to 


p - (h + 1) = 


-(2f - f 2 cos 20)p 


_ * (1 ~ cos 40) , 1 * , 5 x2 

8 (h + 1) J 2 + 16 


4(f + f 2 )cos 20 + f 2 cos 40 
& 16 


Again, me can assume that the difference between the quantity h + 1 and the geo- 
centric distance p can be expanded in powers of the earth's flattening f. This power 
series expansion can be written as 
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(56) 


p - (h + 1) = bj’ f + b 2 ' f 2 

where the terms of order f® and higher have been omitted and the coefficients bj ’ 
and b 2 ’ are to be determined. 

Substituting equation (56) into the left-hand side of equation (55) and then equating 
like powers of f yield the coefficients bj' and b 2 * as follows: 

bj* = — 2 f ~ 1 (57) 

b 2’ = [ioTTTj + 1| (1 ‘ cos 401 < 58 > 


Substituting equations (57) and (58) into equation (56) proves that 


p = (h + 1) + 


( cos 20-1 
2 




- cos 40)) V 


(59) 


Therefore, equation (59) is the transformation formula, in the form of a series expansion 
in powers of the earth's flattening, for the geocentric distance as a function of the geo- 
detic latitude and the altitude. 

Summation of formulas .- The four transformation formulas for the geodetic latitude 
and the altitude as a function of the geocentric latitude and the geocentric distance for 
nonzero altitudes and vice versa, in the forms of series expansions in powers of the 
earth's flattening, are given by equations (36), (46), (53), and (59), respectively. For con- 
venience, these formulas are listed as follows: 




_1 + _J_1 

1 + l)2 + 4(h + l)J 


0 ' = <p + 


( -sin 2<A , J -sin 20 

h + l /|2(h + l) 2 + 




sin 4 <p Vf 
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These power series expansions in f are useful forms for two reasons: (1) the power 
series expansions allow one to readily discern the accuracy of the formulas, and (2) the 
earth’s flattening f is an important physical parameter. 

Transformation Formulas in Forms of Series Expansions 
in Powers of Earth's Eccentricity 

The earth's eccentricity is also a useful and important physical parameter. Equa- 
tion (4) relates the eccentricity e to the flattening f. From equation (4) 

f = 1 - (1 - e2) 1 /2 

Expanding this equation in a binomial series expansion and omitting the terms of 
order e 6 and higher establish that 

f = I ® 2 + 8 ® 4 (60) 

Substituting equation (60) into equations (36), (46), (53), and (59), respectively, and 
omitting the terms of order e® and higher give 


h = (p - 1) + /1 jlCOs 2^\ e 2 + [ ~(1 - cos 20’) / 1 l\ 

\ 4 / (_ 16 \16p 64 ) 


0 ’ = 0 + 


P k + — - - — vjsin 4<^ 

[l6(h + l) 2 I6(h + l)j 


(1 - cos 40’) 


(61) 


(62) 


( 63 ) 




Therefore, equations (61), (62), (63), and (64) are the four transformation formulas, in 
the forms of series expansions in powers of the earth’s eccentricity, for the geodetic 
latitude and the altitude as a function of the geocentric latitude and the geocentric distance 
for nonzero altitudes and vice versa. Note that equation (61) for the geodetic latitude, if 
rearranged, is the same as the Fourier series expansion for the geodetic latitude found 
in reference 9. 


Transformation Formulas in Forms of Fourier Series Expansions 

Fourier series expansions are also useful forms. Equations (36), (46), (53), 
and (59) can be rearranged to represent Fourier series expansions in terms of the geo- 
detic latitude or the geocentric latitude. Therefore, the four transformation formulas 
for the geodetic latitude and the altitude as a function of the geocentric latitude and the 
geocentric distance for nonzero altitudes and vice versa, in the forms of Fourier series 
expansions, are given by the following equations: 




(65) 


( 66 ) 


(67) 


(68) 
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Note that equation (67) for the geocentric latitude, if the altitude h is set equal to 
zero, reduces to the Fourier series expansion for the geocentric latitude found in refer- 
ences 3 to 5. Also, note that equation (68) for the geocentric distance, if the altitude h 
is set equal to zero, reduces to the Fourier series expansion for the geocentric distance 
found in references 3 to 5. Therefore, in this paper, the Fourier series expansions for 
the geocentric latitude and the geocentric distance as a function of the geodetic latitude 
have been generalized to include nonzero altitudes. 

CONCLUDING REMARKS 

Four formulas, for the nonzero altitude transformation from geodetic coordinates 
(geodetic latitude and altitude) to geocentric coordinates (geocentric latitude and geocen- 
tric distance) and vice versa, are derived. The set of four formulas is expressed in each 
of three useful forms: series expansion in powers of the earth’s flattening; series expan- 
sion in powers of the earth's eccentricity; and Fourier series expansion in terms of the 
geodetic latitude or the geocentric latitude. The error incurred due to omitting the third 
order and higher terms in the series expansions is of the order of one part in 3 x 10*, 
which for most environmental and space sciences applications is inconsequential. This 
paper has three significant contributions: (1) the four transformation formulas, which 
are accurate for any altitude, have been derived as series expansions; (2) the formulas 
are listed in a complete and consistent set; and (3) the set of formulas has been cast into 
three different and useful forms. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 6, 1974. 
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